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Abstract 

In this paper, we consider a Partial Interference Alignment and Interference Detection (PIAID) 
design for if -user quasi-static MIMO interference channels with discrete constellation inputs. Each 
transmitter has M antennas and transmits L independent data streams to the desired receiver with 
N receive antennas. We focus on the case where not all K — 1 interfering transmitters can be 
aligned at every receiver. As a result, there will be residual interference at each receiver that cannot 
be aligned. Each receiver detects and cancels the residual interference based on the constellation 
map. However, there is a window of unfavorable interference profile at the receiver for Interference 
Detection (ID). In this paper, we propose a low complexity Partial Interference Alignment scheme 
in which we dynamically select the user set for IA so as to create a favorable interference profile 
for ID at each receiver. We first derive the average symbol error rate (SER) by taking into account 
of the non-Guassian residual interference due to discrete constellation. Using graph theory, we then 
devise a low complexity user set selection algorithm for the PIAID scheme, which minimizes the 
asymptotically tight bound for the average end-to-end SER performance. Moreover, we substantially 
simplify interference detection at the receiver using Semi-Definite Relaxation (SDR) techniques. It 
is shown that the SER performance of the proposed PIAID scheme has significant gain compared 
with various conventional baseline solutions. 
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I. Introduction 

Interference has been a very difficult problem in wireless communications.. For instance, the 
capacity region of two-user Gaussian interference channels has been an open problem for over 30 years 
[1]. Recently, there are some progress made in understanding the interference Q, Q and extensive 
studies have been done regarding the Interference Alignment (IA) Q. For instance, IA is a signal 
processing approach that attempts to simultaneously align the interference on a lower dimension 
subspace at each receiver so that the desired signals can be transmitted on the interference-free 
dimensions. In Q, the authors show that IA (using infinite dimension symbol extension in time or 
frequency selective fading channels) is optimal in Degrees-of-Freedom (DoF) sense. In |5] the authors 
propose a variation of the IA scheme, the ergodic alignment scheme, for K-user time or frequency- 
selective interference channels. In practice, since it is not possible to realize infinite dimension symbol 
extensions, there are a number of works O-lEl and the reference therein that consider IA in the 
spatial domain, without symbol extensions, in the if-user quasi-static MIMO interference channels. 
Furthermore, the authors in |9] investigate the asymptotic performance of these different IA solutions. 
However, for quasi-static (or constant) channels, conventional IA might be infeasible depending on 
the system parameters. For example, it is conjectured in iflOll that conventional IA on quasi-static 
MIMO (M transmit and N receive antennas) interference channels is not feasible to achieve a per 
user DoF greater than ^jnp ■ As a result, we cannot rely on IA to eliminate all interference in quasi- 
static MIMO interference channels especially when K is large. There are other works that consider 
IA over signal scale in |[ITI - lfT4ll . In ifTTll and Ifl2l . signal scale alignment schemes are introduced 
for the many-to-one interference channel and fully connected interference networks, respectively. In 
lfT3l . lfl4ll . the authors propose a lattice alignment scheme for K-user MIMO interference channels. 
However, the scheme requires infinite SNR and serves only as a proof of concept. It is not clear 
whether this approach can be applied at finite SNR. 

Due to the fact that not all the interferers can be aligned at each receiver, there will be residual inter- 
ference at the receiver. In this paper, we assume the receiver has interference detection (ID) capability. 
Specifically, the receiver detects and cancels the residual interference based on the constellation map 
derived from the discrete constellation inputs. However, there is a window of unfavorable interference 
profile for ID at the receiver. For instance, ID at the receiver is more effective when the interference 
is stronger than the desired signal lfl4l . |[T5l . In lfl31l . irT6l . the authors propose an ID scheme for 
quasi-static interference channels based on lattices. However the proposed scheme can only work 
under idealized assumptions such as the symmetric SISO interference channels (where all cross links 
have the same fading coefficients) or a specialized class of 3-user SISO interference channels (where 
the product of fading coefficients is assumed to be rational). 
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In fact, IA and ID are complementary approaches to deal with interference in quasi-static MIMO 
interference channels. The IA approach^ can be used to first eliminate some interference and the ID 
approach can be used to deal with the residual interference at each receiver. While we can potentially 
benefit from the concepts of IA and ID approaches in dealing with interference, there are still some 
key technical challenges to be addressed. 

• Feasibility Issue of IA and Path Loss Effects: Sometimes brute-force IA in quasi-static 
MIMO interference channels (without symbol extensions) might be infeasible depending on 
system parameters. Furthermore, existing literature has completely ignored the effects of path 
loss, which may also be exploited when dealing with interference. For instance, nodes with large 
path loss may not need to be interference-aligned and hence, it is important to jointly consider 
the feasibility issue and the path loss effects. 

• Coupling between IA and ID: While both IA and ID are effective means to mitigate interference, 
their designs are coupled together in an intricate manner. For instance, the performance of ID 
at the receiver depends heavily on the interference profile. Hence, IA can potentially contribute 
to creating a more desirable interference profile for ID by careful selecting a subset of users 
for IA. However, the problem of user selection in IA to optimize the symbol error rate (SER) 
performance is very complicated. First, the optimization space is combinatorial and brute force 
exhaustive search is not viable. Second, even if we can afford the search complexity, obtaining 
the search metric is highly non-trivial because it is also very challenging to analyze the closed 
form average SER under non-Gaussian interferences. 

In this paper, we propose a low complexity Partial Interference Alignment and Interference Detec- 
tion (PIAID) scheme for .fT-user quasi-static MIMO interference channels with discrete constellation 
inputs. We consider QPSK constellations^ at the inputs of the K transmitters and each transmit-receive 
pair may have different path losses. The proposed PIAID scheme dynamically selects the interference 
alignment set at each receiver based on the path loss information to create a favorable interference 
profile at each receiver for ID processing as illustrated in Fig. |2] Interference alignment is applied 
only to the members of the alignment set. We derive the average SER by taking into account the 
non-Guassian residual interference due to discrete constellation. Using graph theory, we transform 
the combinatorial problem into a linear programming (LP) problem and obtain a low complexity user 
set selection algorithm for the PIAID scheme, which minimizes the asymptotically tight bound for the 
average end-to-end SER performance. Furthermore, using Semi-Definite Relaxation (SDR) technique 



1 In the remaining papers, the mentioned IA approach specifically refers to the signal space alignment approach 

Q 

easily. 



2 QPSK constellations are easier to analyze. However, the proposed framework can be extended to QAM constellations 
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ifTTl - EOl . we propose a low complexity ID algorithm at the receiver. The SER performance of 
the proposed PIAID scheme is shown to have significant gain compared with various conventional 
baseline solutions. 

Outline: The rest of this paper is organized as follows. In Section [TTJ we outline the system model 
and the proposed PIAID scheme. In Section|IIIl we discuss the optimization of the IA user selection. In 
Section JVj we derive the average SER by taking into account the non-Guassian residual interference 
and propose a low complexity ID algorithm at the receiver that uses SDR technique. The numerical 
simulation results are illustrated in Section |Vl Finally we conclude with a brief summary in Section 

ED 



A. K-user Quasi-Static MIMO Interference Channels 

We consider i'T-user quasi-static MIMO Gaussian interference channels as illustrated in Fig. [TJ 
Specifically, each M-antenna transmitter, tries to communicate to its corresponding iV-antenna re- 
ceiver. The channel output at the fc-th receive node is described as follows: 



where K = {!,■■■ ,K}, H ki G C NxM is the MIMO complex fading coefficients from the z-th 
transmitter to the A;-th receiver, L^i is the long term path gain from the i-th transmitter to the /c-th 
receiver, and Pj is the average transmit power of the i-th transmitter. In £[)), Xj G C Mxl is the 
complex signal vector transmitted by transmit node i, and z& G C Nxl is the circularly symmetric 
Additive White Gaussian Noise (AWGN) vector at receive node k. We assume all noise terms are 
i.i.d zero mean complex Gaussian with Efz^Zfc)^] = 21 jv- Furthermore, the assumption on channel 
model is given as follows: 

Assumption I (Assumption on Channel Model /T271/): We assume that the long term path gain is 
given by L ki = wcfT. 7 , where d^i is the distance between transmit node i and receive node k, uj is the 
Log-normal shadow fading with a standard deviation a u , and 7 is the path loss exponent. Furthermore, 
we assume that the entries of H^j for all k, i are i.i.d. complex Gaussian random variables given by 
[Hfcj]( n m ) ~ CAA(0, 1) for all k,i,n,m, where [Hfci]( n , m ) denotes the (n, m) th element of H^. ■ 

In this paper, we assume that the i-th transmit node transmits D < min(M, N) independent QPSK 
data streams {x}, - ■ ■ , xf } to the z-th receive node where xf G S = {^(1 ^(1 - j), ^(— 1 + 
j)i 2 ( — ^ ~ 3')}iVd G {1, • • • ,D}. Let (||v^|| = 1, where || • || denotes the Frobenius norm.) 
denote the precoder for the xf symbol. Hence, the transmitted vector at the i-th transmitter is given 



II. System Model 




(1) 



b y x * = Ed v f^- 
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B. Partial Interference Alignment and Interference Detection (PIAID) Scheme 

The proposed PIAID scheme consists of two major components, namely the Partial Interference 
Alignment (PIA) at the transmitters and the Interference Detection (ID) at the receivers as illustrated 
in Fig. |2] 

1) Overview of PIA: PIA is motivated by the feasibility issudjof the MIMO interference alignment 
without symbol extension ifTOll . For instance, it is conjectured in IfTOl that only when K < 
M±K _ x, the K — 1 interfering transmitters can be aligned at every receiver node. As a result, not 
all the K — 1 interfering transmitters can be aligned at every receiver node for large K. Furthermore, 
existing IA schemes do not consider or exploit the effects of different path losses between transmit 
and receive pairs. When the path loss effects are taken into consideration, not all transmitters will 
contribute the same effect at the receiver and hence, there should be different priority in determining 
which nodes should be aligned given the feasibility constraint. Fig. [3] illustrates an example of 4-user 
interference channels where M = N = 2 and D = 1. Using the feasibility condition of MIMO IA 
in ifTOl . only two transmitters can be aligned at each receiver. Combining with the path costs (which 
depend on the path gains and transmit power), it is obvious that transmitters 2 and 3 should be aligned 
at receiver 1 as indicated in the Fig. [3] 

Motivated by the above example, the proposed PIAID scheme dynamically selects a transmitters 
to be aligned at each receiver node based on the path costs. The index of the aligned transmitters at 
each receiver is given by a PIA set with cardinality a. Specifically, the PIA set is defined below. 

Definition 1 (PIA Set): A PIA set is defined as A = {Ak, V/c}, where Ak = {k±, &2, • • • > k a ^ k} 
denotes the index of aligned transmitters at the receiver k for some constant a. ■ 

Only the transmit nodes that belong to Ak will have to align their transmit signals by choosing 
the precoders and equalizers according to the traditional IA requirement! IfTOll . IT221 

(U fe )tH fei Vi = 0,VieA k 

(U fc )tH fefc V fe = diag(Ai,--- ,\ D ),Vk (2) 

where (-)t denotes the Hermitian transpose, diag(Ai, • • • , Ad) denotes a diagonal matrix with diagonal 
entries Ai, • • • , Xd, Ufc = [u^,, • • • , u^] are the N x D decorrelators at receiver k with ||u' fc || = 1, V/, 
and Vj = [vj, ■ ■ ■ , wf] are the M x D precoders at transmitter i with ||v'|| = 1,VZ. Based on the 

3 For constant MIMO interference channels, it is not always possible to completely align all the K — 1 interferers at 
each receiver. 

4 Note that the requirement is equivalent to that used in 1101 . 1221 . i.e., choosing the precoders and equalizers satisfying 
rank((U fc ) t H fefc V fc ) = D,Vk. 
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conjecture in ifTOl . a sufficient condition^ for a feasible PIA set A is given by ifTOl : 

a(A) <mm(^ M p N - 2,K-l) and l(i 6 A k ) = a(A), Vi (3) 

where l(-) is the indicator function, a(„4) = is the cardinality of Ak (i-e-, the number of aligned 
interferers at each receiver). The requirement as per (0 means that each transmitter should be selected 
by a(A) receivers. The IA set A is a design parameter in the proposed PIA scheme, and how to 
choose the IA set is presented in Section [TTTJ 

Remark 1 (Feasibility Condition): Note that the feasibility conditions for ((2]) are still open in the 
literature. Since the feasibility condition and computation of {XJk, Vfc,V£;} are not the focus of the 
paper, we have adopted the results in IfTOl to derive a sufficient condition ((3]). While this condition 
restricts the choice of the feasible set A, it introduces graph structure for the optimization w.r.t. 
A. Furthermore, as shown in Fig. [6j the proposed PIA algorithm with condition Q has similar 
performance as the solution obtained by brute-force exhaustive search (without (O). ■ 

2 ) Overview of ID: In this paper, we focus on the case where not all the K — 1 interferers can 
be aligned at each receivers. As a result, there will be residual interference at every receiver. The ID 
processing at the receiver first estimates the aggregate interference signal by using the constellation 
maps derived from the QPSK inputs. The desired signal is detected after subtracting the estimated 
aggregate interference. For instance, the normalized received signal at the /c-th receiver is given by: 

y/c = V PkLkkH-kk v k X k + ^2/^1 \/ PkL kk^-kk^k x k 



inter-stream interference 

+ J2 -a a .a V p i L kiKkivfxf + V ^PiL ki n ki vfxf +z k . 

non-aligned interferers aligned interferers 

We adopt linear processing at the receiver and the detection process for Z-th data stream at the k-th 
receiver is divided into two stages, namely the aggregate interference detection stage (stage I) and 
the desired signal detection stage (stage II). The two stages are elaborated below: 

• Stage I Processing: Using the /-th column of Ufc in (O, ui, as the decorrelator, the post- 
processed signal of the Z-th stream is given by: 



y{ = (u[^yk = ^/PkTk~kH l k l k x l k + £ yfiV^H 1 *** +(u l k )^ k , (5) 



non-aligned interferers 

5 From 1101 . we know that the total number of equations for the IA requirement {(Ufe)^HfciVj = 0, Vi £ Ak,Vk} is 
N e = J2 k IAI-D 2 = aKD 2 , and the total number of variables is N v = J2 k D(M + N ~ 2D) = KD(M + N- 2D). 
When each transmitter is selected by a receivers, the feasibility condition is simply given by N e < N v 1101 . i.e., a(A) < 
min(M±K - 2, K - l). 
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where H ki = ( u i) H&jV^ is the equivalent channel gain for the d-th data stream of transmitter i at 
receiver k. Note that the inter-stream interference and the interference contributed by users in the 
IA set Ak are completely eliminated due to the PI A requirement in (O. Let Qk = {i '■ PiLki > 
PkLkk^i Ak,i 7^ k} C {1, • • • , K} denotes the set of strong residual interference. Since ID 
at the receiver is more effective when the interference is stronger than the desired signal lfl4T . lfl31 . 
the first stage processing estimates the aggregate strong interference I l k = J2ieQ k -d V r &iLkiHj£xf 
using the following nearest neighbor detection rule. 

Algorithm 1 (Stage I Interference Detection Algorithm): Based on the decorrelator output, yi, 
the detected aggregate strong interference P k is given by: 

(4)* = argrm n bfe- c l> (6) 

where Z l k = { Yl V ' Pi^kiH^sf : sf G S] is the set of possible values the strong interference 
from Qk can take. ■ 
Note that when Qp. = 0, there will be no stage I decoding for the desired data stream x l k . In this 
case, the proposed PIAID scheme reduces to the conventional receiver with one-stage decoding. 
Stage II Processing: The estimated aggregate strong interference (I l k )* is first subtracted from 
the decorrelator output yi as illustrated: 



y 4 k = y{- (4)* = VP k L kk H l k l k x l k + {![)* + ^/PO^H l k ixi+(u l k )^ k , (7) 



ieO k :d 



residual interference 

where O k = {i : PiLki < PkLkk^i Ak,i 7^ k} C {1, ••• , K} denotes the set of weak 
residual interference, and obviously we have Ok [J Qk U -^k U{^} = m turn ' the desired 
signal for receiver k is detected based on y l k using the following algorithm. 
Algorithm 2 (Stage II Signal Detection Algorithm): The Z-th data symbol at the k-th receiver x l k 
is detected based on y l k according to the minimum-distance rule given by: 

& if a( f I >o , if 9/ I >o 

else 



Note that the performance of the ID processing depends heavily on the interference profile, which 
contains the relative power of the residual interference at the receiver. Fig. |4] illustrates the average 
end-to-end SER performance of the ID detector versus the interference power. Observe that there 
is a window of unfavorable interference power for which the performance of the ID is quite poor. 
As such, the user selection of the PIA stage can contribute significantly to the end-to-end SER 
performance of the PIAID scheme. Intuitively, the user set selection A of PIAID should not aim at 
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removing the strongest interference. On the contrary, it should remove the unfavorable interference 
characterized by the ID stage requirement (similar to Fig. [4]). As a result, the PIA and ID processing 
are complementary approaches to combat interference and their designs are tightly coupled together. 

III. Dynamic IA Set Selection in Partial Interference Alignment 

In this section, we shall formulate the dynamic selection of PIA set as a combinatorial optimization 
problem, and derive a low complexity optimal solution using graph theory and linear programming. 

A. Problem Formulation 

Due to heterogeneous path losses and transmit powers, interference links have different contributions 
to the average end-to-end SER of the PLAID scheme. We shall first formulate the dynamic PIA set 
selection problem using a general cost metric C = {cfcj,V/c,i}. In the next section, we shall obtain a 
specialized cost metric related to the end-to-end SER of the PIAID. The PIA set optimization problem 
is summarized below. 

Problem I (MaxPIA Problem): Given a general cost matrix of the interference links C = {cki, Vfe, i}, 
the MaxPIA prolbem is given by: 



where A* is the solution to the MaxPIA problem, and A denotes the collection of all the PIA sets 



B. Optimal Solution of the Dynamic PIA Set Selection 

Optimization problem © is a constrained combinatorial optimization problem, which is difficult in 
general. Solving problem ([9]) using brute force exhaustive search has a high complexity of 0(exp(K)) 
and is not viable in practice. In this section, we shall exploit specific problem structure and visualize 
the optimization problem in (© using graph theory. 

We first review some preliminaries on graph theory from ll23l . ll24l and the reference therein. A 
graph Q is defined by a pair Q = (W,£), where W is a finite set of nodes and £ is a finite set of 
edges. Specifically, the nodes in W are denoted as w\,W2,--- , w\yj\, and an edge in £ connecting 
nodes wi and Wk is denoted as [wi,Wk]. If an edge e = [101,102] G £, then we say that e is incident 
upon wi (and 102). The degree of a node w of W is the number of edges incident upon w. A bipartite 
graph is a graph B = (W,£) such that W can be partitioned into two sets, W\ and W2, and each 
edge in £ has one node in Wi and one node in WV The bipartite graph is usually denoted by 
B = (Wi, W2,£). An example of a bipartite graph is illustrated in Fig. [3] 



MaxPIA: 




(9) 



that satisfies the IA feasibility condition ©. 
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In fact, interference networks can be represented by a bipartite graph B = (1Z,T,£), where TZ = 
{fi, • • • , tk] is the set of the receive nodes, T = {t\, ■ ■ ■ , t^} is the set of the transmit nodes, and 
£ = { [rjfc, tj], VA;, i and k / i} is the set of the edges. A feasible PIA set A is equivalent to a subset of 
the edges £ s with the property that the degree of each receive and transmit node of Q s = (T,1Z,£ S ) 
is a(A), and Q s is called a a-factor of graph Q. 

Example I (Graph Illustration of the PIA set): Suppose the PIA set is given by A = {A\ = 
{2, 3}, A2 = {3, 4}, As = {4,1}, A4 = {1,2}}, the corresponding subset of edges £ s is given 
by £ s = {[ri,t 2 ], [n,t 3 ], [r 2 ,t 3 ], [r 2 ,U], [r3,U], b"3,*i], b"4,*l], b"4,*2]} as illustrated in Fig. [3] ■ 

Let Cki denote the cost of edge [r^, ij] 6 £. Problem [T] is equivalent to finding a a-factor of Q with 
the largest sum of costs. Hence, the MaxPIA problem in © is similar to a matching problem (finding 
a "best" 1 -factor of graph Q) on a bipartite graph, and exploiting this equivalence we shall derive 
a low complexity optimal solution. Let e = {efcj,VA;,i} be a set of variables. If the edge [r^,^] is 
included in the a-factor (i.e. the transmit node i is chosen as one of the aligned interferers at receive 
node k) then e&f = 1, otherwise e^j = 0. As a result, the problem Q] is equivalent to 

argmax £fci Y,k,i c ^ki 
s-t Y,i e ki = a,\/k 



KJ 



(10) 

e k i e {0, l},Vk,i 



where a = min (^±^ - 2, K - l). 

The above problem is a non-convex problem due to the non-convex constraint G {0, 1}. To get 
a low complexity solution, we first relax the constraints eki £ {0, 1} to < eki < 1. As a result, 
(fTOl) becomes a standard LP problem, which can be solved efficiently by the well known simplex 
algorithm [23], |[24l . The following Lemma summarizes the optimality of this relaxation. 

Lemma 1 ( Optimality of the LP Relaxation ): The optimal solution of the LP relaxation problem is 
also the optimal solution of (fTOl) . i.e. et i £ {0, 1} where eL is the optimal solution of the LP relaxed 
problem. ■ 
Proof: Please refer to Appendix [Aj ■ 

IV. SER Analysis and Low Complexity ID Processing 

In this section, we first derive the average SER of the PIAID scheme for a given PIA set A and 
the path gains {L^Pi, VA;, i}. Based on the SER results, we obtain an equivalent cost metric for the 
PIA set selection optimization problem in (fTOl . which is order-optimal w.r.t. the average end-to-end 
SER. Finally, we propose a low complexity ID algorithm using SDR technique. 
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A. SER Analysis of ID with Non-Gaussian Residual Interference 

Unlike standard SER analysis in existing literature ll25l . |[26l . a key challenge of SER analysis in 
the PIAID scheme is that the residual interference in the stage I and stage II ID processing are non- 
Gaussian due to the discrete constellation inputs. As a result, we focus on deriving an asymptotically 
tight SER expression for the interference dominated regime. Theorem Q] summarizes our main results. 



Theorem I (Average SER of the Two Stage PIAID Processing): For a given PIA set A, the SER 
of the Z-fh data stream at the A;-th receiver in the interference limited regime of the PIAID scheme is 
given by 

SER l k (A) 4 E«[SERi(A^)] =E«[E4PK4}m4^4l4>A^}l 

, \ J (11) 

- ^ yl^ieQt P t L ki ^ l^i£O k P k L kk J> 

where % = {H^, Vfc, i}. g(x) = @(f(x)) denotes lim sup a .^ |M < C and lim inf — c 
for some constants C > c > 0. ■ 
Proof: Please refer to the Appendix |B] ■ 

Remark 2: In this paper, the precoders and decorrelators {Vfc,Ufc,V&;} are determined based 
on the IA requirements in ©. Hence, they are only dependent on the channels in the set A for 
which interference is aligned. Since the interference from these channels involved are nulled, the 
remaining interference has a random channel matrix even though it is now projected on the space 
Ufc. Furthermore, the SER in (TTTT) is averaged over realizations of the channels and noise. ■ 

Remark 3 (Interpretation of Theorem^: The result in (fTTb indicates that the SER of the PIAID 
scheme favors either very strong or very weak residual interference. In other words, there is always an 
unfavorable window of residual interferences as illustrated in Fig. [4] The role of PIA is to eliminate 
these unfavorable windows of interferences so that the ID processing is given a more favorable 
interference profile. ■ 

Motivated by Theorem [TJ we set the interference cost metric in the PIA set optimization problem 

as 

f -C ifi = k 

else if PiL ki > P k L kk ,Vfc G K (12) 

□*f M otherwise 
r k ij kk 

where C > is a large constant (a sufficiently large C can be chosen as: C > ^ fc ijLk \c k i\). Based 
on these interference cost metrics, the PIA set selection solutions solved by the LP relaxation is 



June 28, 2011 



DRAFT 



to 



order-optima 



HI w.r.t. 



the following problem: 



A* = arg min V SER^ (A) . 



(13) 



S. Low Complexity ID 

Note that the complexity of decoding algorithm [TJ in the stage I processing is exponential w.r.t. the 
cardinality of the set of strong residual interferences, i.e., \Q k \- Using the SDR technique tfTTl - EOl . 
we shall derive a low complexity ID algorithm, which has polynomial complexity w.r.t. \Q k \. The 
SDR technique has been widely used in multiuser detection l?T71 . lfT8l and MIMO systems |fl9l , EU1 
to derive low complexity suboptimal detectors. It has been shown that the SDR detector can provide 
better performance compared with other suboptimal detectors. 

To utilize the SDR technique, we first simplify equation © as follows: 

(h Q ) T x Q + y/P k L kk H k l k x l k + £ y/iV^Hgxt + (u£)tz fc , (14) 



ieO k :d 



where xg = \x\ x , • • • ,xP,xj 2 ,--- , xf Q ] T G S D Q xl is DQ interference symbols, and Q = \ Q k | is the 



rW it 



hi 



cardinality of Q k . h Q = ly/P^Hf^ y/^L^H^ , y/^L^Hj^ , • • • , y/P~L k ~H : 
G C D ® xl is the channel gain for the interference symbols xg. Furthermore, the real valued form of 
([141 can be expressed as: 

'•>?{"..}' 
^{ n o}. 







-5?{(h Q n -9{(h Q ) T }" 




"K{xc 


2l" 






9{(h e n »{(h e F}_ 








V v 











(15) 



where y R £ l 2xl , H K G K 2x2£ >Q, Xfi £ { ±1 }2DQxi and n/? g M 2xl . Decoding algorithm [Tj in 
stage I processing is equivalent to 



arg min \\y R - H_ R x R \\, 



(16) 



such that the detected aggregate strong interference (I k )* = (hg) T (xg)*, where (xq)* is determined 
from (x. R )* as indicated in (fT51) . The above problem (fT6l ) can be equivalently expressed as 

arg min trace(WS) 



s 

S.t 



diag(S) = I 2 dq+i 



(17) 



ss 



r 



6 It can be observed from the following fact. From Theorem Q] we have A* = 

argmin^ eA © (j2 k ^ Qk T^f + £ Me0fc = ar S min .AeA © ( E*,<** p?XI7 1 ( P * L '" > PkL kk )- 

l{PiL ki > P fc L fefc ) + E fc ie A k T^t 1 ^" < PkLuk))- Hence, the 

( 



argmax^ A 0(^E fc , !e ^^f' 

Efc.i e Q fc +Efc,i e o fe PfcLfct ) is an asymptotically tight bound for SER when P k L kk > PiLfci or PfcL fcfc < PiL fc i 

for all {k, i}, and an order-optimal solution means that it minimizes the asymptotically tight bound for the SER. 
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where (-) T denotes the transpose, 

(H R ) T H R -(H R fy R 



W 



and 



L 1 



(18) 



-(y R ) T U R 

By means of SDR, we relax the constraint S = ss T by S >z (i.e., S G C nxn is positive- 
semidefinite), and problem (fTTT ) degenerates into the following Semi-Definite Program (SDP) that can 
be solved efficiently in 0(Q 3 5 ) time lfT7l-ll20l. e.g., using the interior-point optimization technique 
E3: 



argmin trace(WS) 

s 

s.t diag(S) = I 2 dq+i 



(19) 



{ sto. 

If the optimal value S* of the above problem (fT9l is rank one, then the relaxation is tight, and the 
optimal solution of the (x#)* is given by lfT71 - lfT9l : 



[( x R)*](n,i) = /x(r,n) 



if 



1 else 



> 



Vn G {!,• 



1} 



(20) 



where r is the eigenvector of S* associated with the only one non-zero eigenvalue. 

On the other hand, if S* is not rank one, then we shall approximate (x.r)* based on S*. Specifically, 
there are a few standard techniques to determine (xr)*, e.g., Randomization, Rank-1 approximation 
and Dominant eigenvector approximation lfT71 - |[T9l . 

To further improve the quality of the approximation, we propose a SDR-SID algorithm based on 
the dominant eigenvector approximation as follows (and illustrated in Fig. [5]>. 

Algorithm 3 (SDR-SID Algorithm): 

• Step 0: Set active set A = {x^,--- , x?, xj 2 , ■ ■ ■ ,x^ Q } that contains all the decoding data 
streams, and the cardinality is A = |A| = DQ. 

Repeat 

• Step 1: According to the active set A, solve optimization problem ( fT9l to obtain S*. 

• Step 2: If S* is rank one, determine (x#)* from (l20l and terminate. 

• Step 3: Extract the A dominant eigenvectors of S*, {ri, • • • , r\}. Compute [xj]( n ^ = / x (rj, n) 
from <HB V* G {1, • • • , A} and Vn G {1, • • • , 2A}. 

• Step 4: Compute Sj = [xj, 1] T [xj, 1], Vz G {1, •*" Choose (x^)* = Xj», where i* = 
arg mm ig{1 ... A} trace(VFSj). 

• Step 5: Determine xf, from (x R )*, where {i*,d*} = argmax i d \\/PiL ki H l ^\l(xf G A). 

• Step 6: Set y\ = y\ — ^/Lk^H l ^xf, , delete x% from the active set A and set A = A — 1. 
Until the active set A is empty. ■ 
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Example 2 (Illustration of Algorithm]^: Suppose y\ in (fl4l ) is given by y\ = \/PLrj^H k \x\ + 
\JP~ 2 L^ 2 H k \x\ + n with \y^L^H l ^ 2 \ > \y/PiLkiH l k \\. The details of the implementation of 
Algorithm [3] are given below. 

• Step 0: Set active set A = {xj;,^} and A = 2. 

• Step 1: Suppose S* is not rank one. Go to step 3. 

• Step 3: Extract the 2 dominant eigenvectors of S*, {ri,r-2}, and obtain {xi,X2} from (1201) . 

• Step 4: Suppose trace(WSi) < tmce(WS 2 ), we choose (xr)* = xj.. 

. Step 5: Since \^P 2 L k2 H l k \\ > | \[P^L k ~\ H l k \ \ , we determine x\ from (x R )*, i.e., x\ = [(xjj)*] (2) i) + 
j[(xr)*]( 4i1 ), which by definition xq = [x^xJ,]. 

. Step 6: Set y[ = y\ - P 2 L k2 H l k \x\ = V P\L kl H l k \x\ + n , A = {x\} and A = 1. Repeat 
from step 1 to obtain x\. ■ 

Algorithm [3] is motivated from the intuition that the error probability of decoding symbol xf,Vi £ 
Qfc,Vd is small if its channel gain y/pLfaH 1 ^ is large. Note that the complexity of Algorithm [3] is 
mainly determined by the complexity of solving the optimization problem ( fT9l ) in the step 1. Since 
the complexity for step 1 is only in 0(Q 3,5 ) time iTTTl - EOl . the overall complexity of Algorithm [3] 
is 0(Q 4 5 ). Furthermore, it can be easily generalized to other approximation techniques by simply 
modifying the way to determine Xj in step 3. Finally, using simulation, we illustrate that the average 
end-to-end SER performance of the low complexity SDR-SID algorithm is similar to the performance 
of Algorithm Q] 

V. Simulation Results and Discussions 

In this section, we evaluate the performance of the proposed PIAID scheme via numerical simu- 
lations. In particular, we compare the performance of the proposed schemes against various baseline 
schemes: 

. Proposed Scheme 1: PIAID with Algorithm [1] (PIAID Algl) 

- The PIA set optimization stage tries to align the unfavorable interference links by setting 
the interference cost metric according to (IT2~T ). 

- ID processing is adopted at each of the K receivers after PIA. Specifically, Algorithm Q] 
and |2] are applied at stage I and stage II processing, respectively. 

. Proposed Scheme 2: PIAID with SDR-SID (PIAID AIg3) 

- The PIA set optimization stage tries to align the unfavorable interference links by setting 
the interference cost metric according to (fT2l ). 

- ID processing is adopted at each of the K receivers after PIA. Specifically, Algorithm [3] 
and |2] are applied at stage I and stage II processing, respectively. 
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• Baseline 1: Randomized PIA (Randomized PIA) 

- The PIA set A is chosen randomly from A, i.e., the collection of all the PIA sets that 
satisfies the IA feasibility condition (f3]). 

- ID processing is adopted at each of the K receivers after PIA. Specifically, Algorithm Q] 
and |2] are applied at stage I and stage II processing, respectively. 

• Baseline 2: Iterative interference alignment (Iterative IA) O Algorithm 1] 

- Alternating optimization is utilized to minimize the weighted sum leakage interference. 

- Conventional one-stage decoding is adopted at each of the K receivers by treating all the 
interference as noise. 

. Baseline 3: Maximizing SINR (Max SINR) (6j Algorithm 2] 

- Alternating optimization is utilized to maximize the SINR at the receivers. 

- Conventional one-stage decoding is adopted at each of the K receivers by treating all the 
interference as noise. 

• Baseline 4: Maximizing Sum-Rate (Max Sum-Rate) Q 

- A gradient ascent approach combined with the alternating optimization is utilized to maxi- 
mize the sum-rate of the receivers. 

- Conventional one-stage decoding is adopted at each of the K receivers by treating all the 
interference as noise. 

• Baseline 5: Minimizing Mean Square Error (Min MSE) (H 

- A joint design to minimize the sum of the MSE of the receivers. 

- Conventional one-stage decoding is adopted at each of the K receivers by treating all the 
interference as noise. ■ 

In the simulations, all the transmit and receive nodes are assumed randomly distributed in a 2km x 
lkm rectangular area as shown in Fig. [T] The channel model is given by Assumption [TJ Specifically, 
we set the log-normal shadowing standard deviation as a w = 12 dB and the path loss exponent as 
7 = 6 as per ITU-R recommendation M.1225 |2T1 . Each transmitter delivers a single stream (D = 1) 
of QPSK symbols. The transmit power of each node is assumed to be the same. 

A. Performance w.r.t. Receive E s /Nq 

Fig. [6] and Fig. [7] illustrate the average end-to-end SER performance per data stream versus receive 
E s /No (dB) of the desired link for K = 5 and K = 6 respectively. The number of transmit and receive 
antennas is given by {M = 3,N = 2}. The number of aligned users for the feasible interference 
alignment is a = 3. The average SER performance is evaluated with 10 7 realizations of noise, 
complex fading coefficients {Hfcj,Vfe,i} and path loss {Lfcj, V&:, i}. Observe that the average SER of 
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all the schemes decreases as the receive E s /Nq increases, and there is significant performance gain 
of the proposed schemes compared to all baselines, even for low complexity PIAID with SDR-SID 
(Algorithm [3]). The performance gain is contributed by the user selection of the PIA stage that moves 
the ID processing out of the unfavorable interference profile as shown in Fig. 0] Furthermore, it 
can also be observed that PIAID with SDR-SID (Algorithm [3]) has similar performance as PIAID 
(Algorithm []]). Finally, Fig. [6] shows that the proposed PIA algorithm with condition © has similar 
performance as the solution obtained by brute-force exhaustive search (without (O). 

B. Cumulative Distribution Function ( CDF) of the SER 

Fig. [8] and Fig. [9] illustrate the Cumulative Distribution Function (CDF) of the SER per data stream 
with receive E s /Nq = 25dB for K = 5 and K = 6 respectively, where the randomness of SER is 
induced by {H^Vfe, «'} and {Lfcj,V/c,z}. The number of transmit and receive antennas is given by 
{M = 3, N = 2}. The number of aligned users for the feasible interference alignment is a = 3. The 
CDF performance is evaluated with 10 7 realizations of noise, complex fading coefficients {H^^k, i} 
and path loss {Lki, V/c, i}. It can also be verified that the proposed scheme achieves not only a smaller 
average SER but also a smaller SER percentile compared with the baselines. 

VI. Conclusions 

In this paper, we propose a low complexity and novel Partial Interference Alignment (PIA) and 
Interference Detection (PIAID) scheme for i^-user quasi-static MIMO interference channels with 
general irrational channel coefficients. Based on the path loss information, the proposed PIAID scheme 
dynamically selects the IA interferers at each receiver such that it moves the ID processing out of the 
unfavorable interference profiles. We derive the average SER by taking into account the non-Guassian 
residual interference, and obtain a low complexity user set selection algorithm for the PIAID scheme, 
which minimizes the asymptotically tight bound for the average end-to-end SER performance. The 
SER performance of the proposed PIAID scheme is shown to have significant gain compared with 
various conventional baseline solutions. 

Appendix A 
Proof of Lemma [T| 

By introducing the slack variables s^i, the LP relaxation problem of (ITOb becomes 

( 

argmax efci Y,k,i c ki e ki 
s.t Ei e M = «>VA; 

E fe e fe j = a,Vi (21) 

e fci + s ki = 1) Vk, i 

eu > 0, > 0,\/k,i 
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and it is equivalent to the following matrix form 

K} = 



(22) 



where e = [en,-- - ,ei K ,e 2 i, 
variables, b = [a, • • • , a, 1, 



, &KK-, S\l 
IT 



If , C 



Cll 



arg max e c e 
s.t Ae = b 
e ^ 

• • , six, S21, • • • , Sif.fir] T is formed by the optimization 
c\Ki c 2i, • • • j cxif, 0, • • • , 0] T is formed by the 



2A" if 2 

path costs, and the matrix A is given by 



K 1 









llxK 


OlxK " 












OlxK 


llxK •• 


■ Oi x K 




B 


2 KxK 2 










O2XXE' 2 


Ik* 


Ik 2 




OlxK 

Ik 


OixK " 
I K •• 


■ llx-ff 

■ IX 
















Ix 2 



(23) 



where nxm denotes an n x m matrix of zeros and I n denotes annxn identity matrix. 

Note that the feasible set for the LP ([22]) is given by the polytope 1Z(A) = {Ae = b, e >z 0}. If 
all the vertices of TZ(A) are integers, there exists an optimal solution such that all the optimization 
variables are integers, and hence the LP (l22l will always lead to an integer optimum when solved by 
the well known simplex algorithm. Since < l,Vfc, i, if is an integer, we have S {0, 1}. As 
a result, it is equivalent to proving that all the vertices of the polytope 11(A) are integers. We first 
define Totally Unimodular (TUM) as follows. 

Definition 2 (Definition of Totally Unimodular): An integer matrix Z is totally unimodular if the 
determinant of each square submatrix of Z is equal to 0,1, or -1. ■ 

It has been shown in 1231 that if A is TUM, then all the vertices of H(A) are integers. Therefore, 
the proof reduces to proving that A is TUM. Note that the matrix B satisfies the following conditions: 

• B is a {0, 1} matrix with no more than two nonzero elements in each column. 

• Each column contains two nonzero elements that have the same sign, where one element is in 
a row contained in the subset Qi = {1, ■ • ■ , K} and the other element is in a row contained in 
the subset U 2 = {K + 1, ■ ■ ■ , 2K}. 

Therefore, B is TUM [24]. It is easy to verify that the matrix C = [B T ,I] T is also TUM. Let D be 
a square, nonsingular submatrix of A. The rows of D can be permuted so that it can be written as 



D 



E 





F 


I 



(24) 



where E is a square submatrix of C, and possibly with its rows permuted. Therefore, det(D) 
det(B) = ±1, which completes the proof. 
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Appendix B 
Proof of Theorem Q] 

A. Upper Bound of the Average SER Performance 

In this subsection, we shall obtain an upper bound of the average SER SER Z fe („4). Specifically, we 
shall first obtain an upper bound of the SER under a given channel realization, i.e., an upper bound 
of SER l k {A,n). 

1) Upper bound of SER k {A,%): Let % l k = {H l ki ,\/i,d} denotes the set of equivalent channel 
gains after applying the equalizer u l k when decoding the desired signal x k at the receiver k, 5 l k = 
l\ — il denotes the difference between the real and estimated aggregate strong interference, and 
C l k = Siee> fc n VPiLkiHj^xf denotes the aggregate weak residual interference. Specifically, we start 
to consider the two stage decoding separately. 

• Stage II decoding: In stage II decoding, y\ in ([7]) is given by: 

Vk = Vk - h = \fPj^kH[ l k x l k + 5{ + Cl + (ui)tz fc , (25) 
Given % l k and suppose 5ft{4l = the error rate of decoding 5ft{4) is given by E6l : 
Pr * = ^,I l k Jl,Hl} = Pr{K{ 7K L} < 0|»{4> = ^,l[JlK 

= q((^+ })\vm^H» k \) = q ivro&i + + cm ) , 

(26) 

where Q(x) = -4- f°° e" x2 / 2 dx denotes the Q-function, and e~ je = ^SSfjL Similarly, 
we have 

Pr{U{x{} * -^\M{x{} = -4 ,I l k J l k ,H{} = Q |vTO^4*J - + £)}) , 

Pr{3{4} * f |3{4> = f , IIHK) = Q (# \Vm^H» k \ + »{e-^(4 + Q}) , 

Pr{Z{x l k } + = -f J l k J l k ,n{} = Q (f IVTOSA&I -»{e-^(^i + a)}) • 

(27) 

Then the error rate of decoding x' fc is given by: 

p r {4^ 414,4,4,74} 

< Pr{»{4> / M{x{}\M{x l k },I l k ,P k ,H l k } +Pr{Qf{4) / 3{4}I»K}>4>4>^} ( 28 ) 
<2Q(^|v / TO^^|-|^|-|C[|). 

• Stage I decoding: In the stage I decoding, y l k in ([5]) is given by: 

y{ = P k + ^P k L kk H l k l k x l k + C l k + (4)t Zfc . (29) 

Note that when the detected aggregate strong interference I l k / il, an upper bound of the 
probability that the estimated aggregate strong interference is P k is given by 

Pr{i l k \x l k ,i l k ,n l k } < Px{P k / 414,4,74} < 2Q (i4i/2 - yp k L kk H l k l k \ - icii) , eo) 
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where Pr{j|, 7^ 4l x L> Hi} is the conditional probability that I k ^ I l k when I k is transmitted. 
Therefore, we can follow the same way as (l26l ) and (l28l ). which obtain the conditional probability 
that xi / when x k is transmitted. When the detected aggregate strong interference I l k = I l k , 
obviously we have following expression at interference limited regime 



Pr{/£|4,4X} = Pr{4l4>4X) < 2Q (-IV^^^ifcl - ICil) « 2. (31) 

As a result, the probability that the estimated aggregate strong interference is I l k in the stage I 
decoding is given by 

From the upper bounJei^skofe^) Ind^^lie "s^fV^s^tllf&Llbding, respectiv ( e$ 
an upper bound of the SER under a given channel realization SERj.(^4, is given by 

= E4&S Pr { x fe} (E/'ezj, Pr {4) Pr {^L 4l4' 4> ?4j) 

= E4&S Pr { x fe} (Ei'ezjj, E/'ex* Pr {4} Pr {^fcl x i'4'^L} Pr {^fc 4l4'4>^i'^i}) 

< ]ij E*» fc6 S E/< e z< Ej* e x£ (l^ll/ 2 - W p kLkkH l k l k \ - ■ 

Q(f iv^^^l-Ki-icil). 

(33) 

2 J £//?per fotmd of SER l k (A): From (|33]>, the average SER is given by 
SER^)^ [SER} k (A,H)] 

< ik\ E n[j:P k eX' k J2iiexi 1(1^1/2 - IvTO^I - IC^I < 0)- 

i(%vmrkH i k i k \-\8 i k \-\ci\<o)_ 

= A Eijezj E/< e z< Pr |VTO^^ fe | - K£l < 1^1 < 2(|vTC^ fe l - |#|)}. 
Therefore, from the above equation (1341) . we shall discuss the following two cases given by: 
successful stage I decoding I l k = I l k (i.e., 5 l k = 0) and unsuccessful stage I decoding I l k ^ I l k 
(i.e., 6[ ± 0). 

• Successful stage I decoding P k = I l k (5 l k = 0): Note that H l ki = (ujjtHfcjvf, and the aggregate 
weak residual interference C} k is given by 

C[ = Eieo fc ,dV ' Pi L kiH l ki xf = Y,ieO k ,d V^fc^u^H^vfzf 

= EieO* V ^fci Em,n [( U L)1 (l,n) [ H to] (n,m) Erf[ V i ^j] (m,l) ■ 

Since {u^v^Vi G Ofc,Vd} are determined by the channel gains from the aligned links given 
by {Hfcj, Vi G ^4^, V/c} ifTUl . 1221 . they are independent of the random variables {H^, Vi G Ofe}- 
From Assumption [T] [Hfcj]( nm ) ~ CAf(0, 1), and hence given {u z fc , vfxf, G Ofc,Vd}, C[ is a 



(34) 
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r d x d 



x - l(m,l) | 



complex Gaussian variable with zero mean and variance 

a 1 = Eieo k p i L kiY, m ,n\[( u k)%,n)\ 2 (| Ed[ v t x f}( m ,i)\) 2 

= EieO k P i L ki E m ( I Ed[v?af]( TO ,l) |) 2 < D Z ie o k P i L ki E m (Ed I t v 

= D 2 Ei£O k P iLkii 

(36) 

because of ||u^,|| = ||vfxf|| = 1. From ( l k is a complex Gaussian variable with zero mean and 
variance a 2 , and the probability density function of the random variable Z = \-jfn-\ 2 is given 
by fz{z) = (z+i) 2 • Therefore, when 5 l k = the conditional probability is given by: 

Pr {^IVTOMJ - |Cil < 0\n{, vfxf,Vi G O fc ,Vd} = Pr{^ < Z} 

= V(^P + 1) < V( 4D ^;I>L t< + 1) = e (E fe0fc fife) • ^ 

In turn, the probability of the event ^\\J P^L^K k l k \ < \<^ l k \ is given by 

Pr{f l^ra^l - Kil < 0} < (£ 4£ o fc fife) • 08) 

Unsuccessful stage I decoding i| / 4 / 0): Given P k = EieQ h ;d V^LkiHlfxf, where 
xf G 5, is given by: 

4= 4-4 = I^Q fc ;d^^W t H«v?(x--af) (39) 

= EieQ* V^wEm l n[( u *) t ](l,n)t H «](n I m) Ed[vfytf](m,l). 

where rpf = xf - xf G {0, ±\/2 ± jV2, ±V2, ±jV2}. Given {uj^vfy^Vi G Q k ,W}, 6 l k is a 
complex Gaussian variable with zero mean and variance given by 

° 2 = E l eQ k P i L ^E m ,n\[( U i)%,n)\ 2 \Edi V t^t](rn,l)\ 2 

= EteQ k P i L kiEm I Edi V f^f](m,l)\ 2 > P i*L ki , Ei^Q k Era I EdWi^i ](m,l) | 2 , 

where i* = argmin i6S PiL^, and H = {i : ^ d |^| / 0}, since ||u^.|| = 1 and ||vfvf|| = \^f\- 
Therefore, when 5{ ^ 0, the conditional probability is given by 

Prj^lVTOMfel " Kil < |^| < 2(|v^X^^| - |C£|)|u< ,vf</f,Vi G Q fc ,Vd} 
< Pr {0 < |d* | < ^VTO^Hjgfcl} = Pr {0 < |J^| 2 < ^fe)} 
< 



4P fc L tt /q 2 ^ 4P fc L fcfc 
4P fc L fcfc /<7 2 +l 



In turn, the unconditional probability is given by: 

Pr{^|WTO*ft| " ICHI < 141 < 2(|VTO^&| " |Cil)} 



(41) 



(42) 



where C = E u ! )V ^d 



,<*„/,<*! 2 



4/(E ieSt E m I EdW^j 

Finally, substitute ([38]> and (@2]> into ([34]), we have 



is a positive constant. 



SERU-4) < (E ieSfe %fe + E^o k fife ) • («) 
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B. Lower Bound of the Average SER Performance 



In this subsection, we shall obtain a lower bound of the average SER SER fc („4). Specifically, we 



shall prove SER fe (^) > e(#£*i.),Vi G O k and SER fc (^) > 0( 



PiL ki h 



Vi G Qfc, respectively. 



Proof of SER fc („4) > e Vt G O fc : Suppose receiver k has perfect knowledge of all 

the data streams except xf,\/i G O k: \/d. After cancelling the known data streams in stage II 
decoding, y l k in (0 is given by: 



—/ 
Vk 



^P k L kk H l k l k x l k + s/P^iHlixf + (ui)tz fc 
Suppose 5ft{4) = the error rate of decoding 3i{4} is given by (cf. (1261)): 
Pr{K{4}^»{4}|K{xi} = ^} 



(44) 



(45) 



where e je = | ^^^ | . Note that ^{e~ j9 H^xf} G W(0, i<7 2 ), and hence the average error 
rate of decoding $t{x l k } is given by: 

Pr{3?{4) / K{4}|3?{4} = ^}] = Pr {^Iv^y^lT&l + R{ e -'»v^#K} < 
Pr {| x /iTO#» fe | < >/TC«} = i(l- y/ i + P iLk )/P k L kk ) 



E 



PL kI , r> , /-;/.:; 







Pfc 

(46) 



where ft G AA(0, <t 2 ) and |i?[ z J G Rayleigh(cr 2 ). Similarly, we have 



E 



Pr{K{4} / 5?{4}|3?{4) = 



E[Pr{9f{4}^9f{4}|9f{4} = ^} 
Pr{^|v^X^^« | +%{e-i e Vm^H l k d i xf} < 0} = 0(#fe), 



(47) 



E 



Pr{9f{4}^9f{4}|9f{4} = -^} 



Pr{^|/fffi^ fc | - ^{e-^VPlLTiH^xf} < 0} = 0( 



P k L k 



Therefore, given that receiver k has perfect knowledge of all the data streams except xf,\/i G 



O k ,\/d, the average SER of decoding x l k is given by 



SER" 



E 



E Pr{4}Pr{4Ml4} 



(48) 



kk 



and obviously we have SER fe („4) > SER^ = 



PiL k 
P h L k 



, Vi G 0£ 



• Proof of SER^A) > 0( p fc ^ fc ), Vi G Q^: Using similar arguments above and suppose receiver 
k has perfect knowledge of all the data streams except xf,Vi G Q k yd. After cancelling the 
known data streams in the stage I decoding, y{ in (f5]) is given by: 



^Pj^H^xi + y/iV^Hgx* + (4)t Zfc . 



(49) 
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Given K{x' fc } = = the error rate of decoding 3t{x' fc } is given by 

Pr{»{4}^»{4}|»{4},»{^}} 
> Pr{»{^}^»{atf}|»{4},»{x?}^^ 



(50) 



where e je = /^i .^ = r^Piir^ = - ip. The probability density function of 

WPkLkkH 1 ' WPiL ki Hl d : -re J J 



27T. As a result, the average error rate of decoding ^{x l k } is given by 



9 I 

■jffi | is given by /z(z) = mrrp » and $ follows the uniform distribution between and 



E[Pr{»{4}^»{4}|»{4},»{xf}}' 
> Pr {2^P)IV^^I < \VKLTkH k l k \ < cos(0)|v?TOr&|} (51) 

_ fV 4 (co S 2 W-l/(4cos 2 W))P fc L fcfc /P,L fc , , q _ flf flL tM 

(l+P fe L fcfc /(4P l L fel co S 2 ( 1 ?)))(l+cos 2 ( 1 ?)P fc L fcfc /P t L fel ) ^ P,L fcl J" 

Furthermore, given that receiver k has perfect knowledge of all the data streams except xf, \/i G 
Qfc , Vd, the average error rate of decoding x\ is given by 

SM^E[Pr {4^414,*?}] >©(%fef). (52) 

Obviously we have SER^(^) > SER? > ©(%jjf ), Vi G Q fc . 
Finally, from the results of upper and lower bound, we can conclude that 

SErL(-4) ^ E n [SER{(A,H)] 4 ^ + ^ gfe) . (53) 
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Fig. 1. Quasi-static K-usei MIMO complex Gaussian interference channels. Each M-antenna transmitter tries to transmit 
D independent data streams to its corresponding TV-antenna receiver. All the nodes are randomly distributed in the a x b 
rectangular area. 
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Fig. 2. System level block diagram of the proposed PIAID Scheme. 



4-users Interference Channels Aligned Transmitters ai Each Receiver 

(a Bipartite Graph) (a 2-factor) 




O O Path cost: 2 O O Path crol: 1 

Fig. 3. An example of 4-user interference channels where M — N = 2 and D — 1. Using the feasibility condition of 
MIMO IA in 1101 . only two transmitters can be aligned at each receiver. Specifically, the 4-user interference channels are 
represented by a bipartite graph B = (TZ, T,£), where TZ = {n, • • ■ , r,i} is the set of the receive nodes, T = {ti, • • • , t^} 
is the set of the transmit nodes, and £ = {[r^, ti], Vfc, i and k ^ i} is the set of the edges. Furthermore, the PIA set is 
given by A = {Ai = {2, 3}, A2 = {3, 4}, ^3 = {4, 1}, Aa = {1, 2}} and the corresponding subset of edges is given by 
£ s = {[ri,t 2 ], [n.,t 3 ], [r 2 ,t 3 ], [r 2 ,U], [r 3 ,U], [r 3 ,ti], [r4.i1], [taM]}- 
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Fig. 4. Average end-to-end SER performance of the ID detector versus the interference power at the interference limited 
regime. Specifically, the input to the ID detector is given by y — H\x\ + \fP%HiXi, where x\ £ S is the desired signal, 
X2 £ 5 is the interference, Hi,H2 ~ CA/"(0, 1), and the interference power is P%. The average end-to-end SER is given 
by SER(P 2 ) = J2 X1 Pr{*i =/= xi\xi,P 2 }. 
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Fig. 5. Illustration of Algorithm [3] (SDR-SID Algorithm). 
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Fig. 6. Average end-to-end SER performance versus the receive E 3 /No (dB). The setup is given by K — 5 (number of 
users), {M — 3, TV = 2} (number of transmit and receive antennas), D = 1 (number of data stream), and a = 3 (number 
of aligned users for feasible interference alignment). 
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Fig. 7. Average end-to-end SER performance versus the receive E a /No (dB). The setup is given by K — 6 (number of 
users), {M — 3, TV = 2} (number of transmit and receive antennas), D = 1 (number of data stream), and a = 3 (number 
of aligned users for feasible interference alignment). 
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Fig. 8. Cumulative Distribution Function (CDF) of the SER per data stream with receive E 3 /No = 25dB. The setup is 
given by K — 5 (number of users), {M = 3, N = 2} (number of transmit and receive antennas), D — 1 (number of data 
stream), and a — 3 (number of aligned users for feasible interference alignment). 
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Fig. 9. Cumulative Distribution Function (CDF) of the SER per data stream with receive E s /No = 25dB. The setup is 
given by K — 6 (number of users), {M = 3, N = 2} (number of transmit and receive antennas), D — 1 (number of data 
stream), and a — 3 (number of aligned users for feasible interference alignment). 
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